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Abstract 

This work addresses a study of coherent states for a physical system governed by a 
Hamiltonian operator describing the motion, in two dimensional space, of spinless elec- 
trons subjected to a perpendicular magnetic field B coupled with a harmonic potential 
^MojqT 2 . The underlying su(l, 1) Lie algebra and Barut-Girardello coherent states 
are constructed and discussed. Further, the Berezin - Klauder - Toeplitz quantiza- 
tion, also known as coherent state (or anti-Wick) quantization, is also performed. The 
thermodynamics of such a quantum electron gas system is elaborated and analyzed. 
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1 Introduction 



The system of charged quantum particle interacting with a constant magnetic field is un- 
doubtedly one of the most investigated problems in quantum mechanics, mainly inspired by 
condensed matter physics and quantum optics. A family of coherent states (CS) adapted 
to such a system was first proposed in pQ. In [2], the behavior of the transverse motion of 
electrons in an external uniform magnetic field B was considered. A complete set of CS wave 
packets was constructed. These states are non spreading packets of minimum uncertainty 
that follow the classical motion. They are the eigenstates of two non-Hermitian operators 
that annihilate the zero- angular-momentum ground state. The CS basis was used for the 
calculation of the partition function. The Landau diamagnetism and the de Haas- van Alphen 
oscillations were contained in this state setting. 

Afterwards, one can also mention some alternative constructions proposed in [3] and [I]. 
In metal and other dense electronic systems, the electrons occupy many Landau levels. Fur- 
thermore, the kinetic energy levels of electrons in two-dimensional gas correspond to Landau 
levels. In particular in [5], the generalized Gazeau-Klauder CS were extended to systems 
with more than one degree of freedom. There, three different types of these generalized 
CS with symmetries of harmonic oscillator, su(2) and su(l, 1) for the Hamiltonian corre- 
sponding to two-dimensional spinless electrons confined by an isotropic harmonic potential 
in the presence of a constant magnetic field were considered. CS also play an important role 
in non-equilibrium statistical physics that describes the evolution towards thermodynamic 
equilibrium for quantum systems with equidistant energy spectra set in thermostat [BJ. 

Besides, the CS were investigated to obtain Landau diamagnetism for a free electron 
gas [TJ. In j8] generalized Klauder-Perelomov [6J and Gazeau-Klauder j9] CS of Landau 
levels were constructed using two different representations for the Lie algebra /i 4 . In pO] , 
the Landau levels were reorganized into two different hidden symmetries, namely su(2) and 
su(l, 1). The representation of su(l, 1) by the Landau levels then led to the construction of 
the Barut-Girardello CS (BGCS) [11]. Moreover, Dehghani et al [12] developed Glauber two- 
variable CS in different representations based on the action of unitary displacement operator 
corresponding to Weyl-Heisenberg algebra on the Klauder-Perelomov CS of su(l, 1) and 
su(2) algebras, for Landau levels, minimizing the Heisenberg uncertainty relation. More 
recently, Bergeron et al [13] investigated the consistency of CS quantization, (also named 
Berezin-Klauder-Toeplitz or anti-Wick quantization), and led to the conclusion that the 
predictions issued from this type of quantization and those from the canonical quantization 
are compatible on a physical level, for non-relativistic systems, even if these two quantization 
techniques are not mathematically equivalent. 

This work aims at considering Landau levels generated by a physical system engendered 
by a Hamiltonian operator describing the motion, in two dimensional space, of spinless 
electrons subjected to a perpendicular magnetic field B coupled with a harmonic potential 
^MuqT 2 . The underlying su(l, 1) Lie algebra and BGCS are constructed and discussed. 
Further, the Berezin - Klauder - Toeplitz quantization is performed through this paper. The 
thermodynamics of such a quantum electron gas system is elaborated and analyzed. 

The paper is organized as follows. In section 2, we start with the study of a spinless 
electron gas on the (x, y)-space in a magnetic field B with an isotropic harmonic potential by 
using step and orbit-center coordinate operators. In section 3, the su(l, 1) representation in 
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the Hamiltonian quantum states is studied. Furthermore, follows a discussion on the BGCS 
of the Lie algebra su(l, 1) defined on the Hilbert subspaces, the mean values of tS77(l,l) 
group generators, the probability density and the time dependence of the BGCS. Some 
quantum mechanical features deduced from the quantization of a complex plane using these 
states are investigated in section 4. The results issued from this procedure are exploited to 
derive quadrature operator mean values and some quantum optical properties. In the section 
5, main statistical properties for the quantum gas in the thermodynamic equilibrium with 
a reservoir (thermostat) at temperature T are computed and analyzed. The last section is 
devoted to concluding remarks. 

2 Hamiltonian operator of an electron in a uniform 
magnetic field with harmonic potential 

Consider a system of spinless electrons (M, e) living on the (x, y)-space in a magnetic field B 
along the ^-direction. Recall that the eigenstates and eigenvalues of such a system were inves- 
tigated for the first time by Landau [7j. When a harmonic confining potential is introduced 
and the Coulomb interactions are neglected, this system is described by the Fock-Darwin 
Hamiltonian [HI [151 
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where p is the canonical momentum and A is the vector potential. We study the problem 
by considering the transverse motion of the electrons on the (x, ?/)-space. 
In the symmetric gauge 
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the classical Hamiltonian "H in (pQ) becomes 



n(p,r)=n 
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The canonical quantization of this system is obtained by introducing the coordinate and 
momentum operators, denoted here by Ri,Pi, respectively, which satisfy 



(4) 



with Ri = X, Y , Pi = P x ,P y , i = 1,2. Denote the Hamiltonian operator associated with 
([3]) by H and solve the eigenvalue problem 



by performing the following change of variables: 



(5) 



Z = X + iY, P z = -(P x -iP y ) 



(6) 
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and considering the set of step operators denned by [2] 



Mtt « MQ, p- 

7t + = 2P Z + i— -Z, tt_ = 2P Z - i—-Z, (7) 



where the motion along the z-axis is free, with Cl = yuj^ + 4cUq and the cyclotron frequency 
w c = |£. These operators satisfy the commutation relation 

[7T_,7T + ] = 2MQh. (8) 

Next, consider the orbit-center coordinate operators given by 

and 

carrying out the following commutation rules 

[X 3t X±]=±hX ± , [X + ,X_] = 2l 2 , [X ± ,i\ = 0, [* 3 ,1] = 0, X 3 = L Z , (11) 



where I := y is taken as the classical radius of the ground-state's Landau orbit for the 
frequency Q. Also, we have 

[7T 3 ,7r ± ] = ±hn ± , [vr±,l] = 0, [tt 3 ,1] = 0, tt 3 = L z . (12) 

Physically, 7r 3 is a constant of the motion of a free particle, as defined in [16], along 
the z— axis; X± step only the angular momentum m and not the energy [2]. Besides, the 
operators X± and tt± commute through the following relations 

[X±,n ± ]=0, [^ ± ,tt t ] = 0. (13) 

Finally, the Hamiltonian fll]) can be rewritten as 
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n+n - r i + ^) + ^-(i-^\xx + hn 
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(14) 



with the eigenstates \n,m) determined by two quantum numbers: n, associated to the en- 
ergy, and m, to the z— projection of the angular momentum, solving the time-independent 
Schrodinger equation H\n,m) = S n)m \n,m) with the corresponding eigenvalues 

£ n>m = Ml (^n + - |(fi - u c )m. (15) 

3 Representation of su(l, 1) algebra by the quantum 
Hamiltonian states 

3.1 Hilbert space representation 

Rewrite now the step and orbit-center coordinate operators n±, X± with the help of dimen- 
sionless variables as follows: 



x + - — 



i . a 2 - 



1 - 2 2 . 

4 ft 



(17) 



where I is defined as previously, so that they satisfy the following canonical commutation 
relations 



7T_,7r_, 



[*+,*-]. 



Introduce the Lie algebra su(l, 1) [T7] corresponding to the SU(1, 1) group, spanned by 
the three group generators {/Ci, /C 2 , /C 3 }: 



[/Ci,/C 2 ] = -«X 3) [/C 2 ,/C 3 ] = i/C x , [/C 3 ,/Ci] = zX 2 . 



(19) 



Use the raising and lowering generators /C± = K\ ± zX 2 to form the second-order differential 
operators given by 



/C+ := 7T+AL, /C_ := it-X + , /C 3 := — (7r + 7r_ + X+X^ 
which satisfy the commutation relations of the Lie algebra su(l, 1) : 

[/C+,/C_] = -2/C 3 , [/C 3 ,/C ± ] = ±/C±. 



(20) 



(21) 



Then consider the Hilbert space = span{\n, 77T.)} n >o,o<m<n = ©m=o^m, where fj m is 
the subspace defined by f) m := span{|n, m)} n > m > . These operators act as follows: 



7r+|n — 1, m — 1) = A/ra|n, m), 7r_|n, m) = ~~ 1? m — 1); 



(22) 



<£+|n, m) = yjri — m)\n, m + 1), X- \n, m) = y/ (n — m+ l)\n, m — 1), 



(23) 



/C+|n — 1, m) = \/n(n — m)\n, m) 
/C_|n,m) = \J n(n — mj|n — l,m) 

/C 3 |n, m) = — (2n — m + l)|n, m). 



|n, m) 



r(m + r 



r(n — m + l)r(n + 1 



/CI 



\m, m). 



(24) 



The lowest state \m, m) is annihilated by the lowering operator JC-, i.e. 

£_|m,m) = (25) 
while the arbitrary states |n,m) are constructed by means of the first equation of (JH 



(26) 
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3.2 CS for the su(l, 1) algebra 

Introducing the BGCS [H] in the Hilbert subspace S) m as eigenstates of the lowering gener- 
ator /C_ of the Lie algebra su(l, 1): 

KJ[z) m = z\z) m (27) 

where z is an arbitrary complex variable with z = pe 1 ^, O<p<oo,O<0<27r, we can 
represent the eigenstates \z) m as the superposition of the complete orthonormal basis \n,m) 

Of $)m ■ 



\z) m =22(n,m\z)m\n,m). (28) 

n=m 

Using (j27p . (1281) . the third equation in (I24p and the orthonormality relation 

+oo 

^2\n,m)(n,m\=I Sim (29) 



we get 



^n{n — m) 
which can be recursively transformed into 



n,m\z) m = ==(n-l,m\z) m (30) 



(n,m\z)m = \l=r f Z r\ z ™ m (m,m\z) m . (31) 

I [n — m + 1JI [n + 1J 

Coming back to the relation ( 1271) . the following expansion of the state \z) m in the Hilbert 
subspace $) m : 



+oo 



'M E zn ' m \l f? 7nr M) K m > ( 32 ) 

>' 1 (n — m + 1)1 (n + 1) 



ii=ni 



holds. The normalization factor |(m, m|^) OT | can be obtained by normalizing to unity the 
states \z) m and using the relation 



+°° x n i 

V — TW/ r = —lJ2x) (33) 



„ n!T(n + z/ + 1) x 

n=0 v ' 

where I u (2x) is the modified Bessel function of order v. It follows 



l(m>m|z)ml = ^ / m (2N)r(m + i) (34) 



and finally, the eigenstates \z) m become 



\z) m = — . > — In, m) (35) 
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where 7 m (2|z|) is the modified Bessel function of the first kind given by [18j: 



+00 



I ^ 1 2n+m 

u 2 N)=E r( „ +1 ; ( „ +m+1) - 

It is demonstrated (see [H]) that these states satisfy the following resolution of the identity: 

\z) mm (z\dg(z) = Jjj m (37) 

which, in this case, holds on i} m where dg(z) is an appropriate measure. Indeed, performing 
the integrals over the whole complex plane, where z = pe^, p £ [0, 00), <fi £ [0, 2ir) and taking 

dg(z) = -I m {2\z\)K m (2\z\)d 2 z, d 2 z = d{Rez)d{Imz), (38) 

71 

the relation (129]) defined on Sj m leads to the resolution of the identity (137)) by using the 
following integral relation: 

POO 

4 / p 2n ~ m+1 K m (2p)dp = T(n-m + l)T(n + 1) (39) 



where K m (2p) given by 



A m (2p) = I / - (2P ' " 'f P) <*» 
2 sin m7r 



is the modified Bessel function of the second kind [18]. The relation ( 1391) is provided by the 
following integral [T9] : 



[ite(//+l±i/)>0,itea>0]. (41) 

The states |z) m form an overcomplete basis for any allowed value of m. Indeed, from (1321 ). 
setting z = pe 1 ^, where p £ [0, 00), (j) £ [0, 27r), it follows that 



\z) mm {z\dg{z) 

= 4 £ r ^f' -t T P 2 — +1 dpln, m> (n, m| (42) 

n=m ' v ' 

and using (l3~9~j) . we get on the Hilbert subspace $y m : 

p +00 

/ \z) mm (z\dg(z) = V |n,m)(n,m| = i^ m . (43) 
./:■ 

n=m 

The BGCS (1351) are not orthogonal. Indeed, given two vectors \z) m and |z') m ' on the Hilbert 
subspaces S) m and $) m > (m 7^ m'), respectively, we get the following expression 

(\z'z\\ ml2 ImhVz 7 ^) 

m (z'\z) m = ^ - V ) =. (44) 

V zfz ) ^I m (2\z>\)I m (2\z\) 
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The overcompleteness of the BGCS \z) m on S) m suggests to discuss their relation with 
the reproducing kernels [20] . 

Define the quantity K.(z,z') := m {z'\z) m on S) m . Using the facts 



n\Z \Z/ m 



\zz> J y /I m (2\z\)I m (2\z'\) ±1 V 



(2n-m)/2 



T(n-m+ l)r(n+ 1) 



(45) 



/lm/2 



m/2 



Am/2 



and setting n = m + i/, it follows that 



(46) 



jn\Z \z^ m 



zz'\\ m/2 I m (2V^) 



zz' J y/I m (2\z'\)I m (2\z\) 
/C(z, z') is well a reproducing kernel. Indeed, 

Proposition 3.1 The following properties 
(i) 



!C(z',z). 



JC(z,z') = K{z',z) : hermiticity; 



(47) 



(48) 



(ii) 



JC(z, z) > : positivity; 



(49) 



dg(z")JC(z, z")JC(z", z') = fC(z, z') : idempotence 



are satisfied by the function K on Sj m . 



(50) 



Proof, (i) and (ii) follow from ( )47|) and ( l44l) . respectively. The left-hand side of ( |50l) is 
written as: 



dg(z")K,(z, z")K(z" , z') = I dg(z") 



\z"z\\z'z"\ 



m/2 



X 



I m (2VWz~) I m (2Vz^) 



^_I m (2\z"\)I m {2\z\)^I m (2\z>\)I m (2\z"\) 
-K m {2\z"\)I m {2\z"\)d 2 z'^ 



7T 



X 



X 



Z'z\\ m/2 (\Z"\\Z"\Y 12 I m (2^) Im QVWF) 



z'z 



z z 



m/2 



I m (2\z"\)y/I m {2\z\)y/I m (2\z'\) 



z'z J y/I m (2\z'\)I m (2\z\) 



X 



x f d 2 z"-K m {2\z"\)I m (2\f¥^j I m (2V¥7^j 

\z'z\ \ m/2 J -( 2 ^) r„..2 



z'z J ^I m (2\z'\)I m (2\z\) Jc 



d 2 z"-K m (2\z"\)I m (2\z"\) 



7i 



= JC{z,z') (51) 
where the relation 

dg(z) = [ d 2 z-K m (2\z\)I m (2\z\) = 1 (52) 
Jc 71 

is used. □ 
From (1431) . for any \^/) G Sj m , we have 



|tf> = / d e (z)V(z)\z) m (53) 
Jc 

where ^(z) := m (z|^). The following reproducing property 

y(z) = I d e {z')K{z, z')^{z') (54) 
Jc 

is also verified. 

The Hilbert space $j m := span{\n, m)} n > m , m > can be represented as the Hilbert space 
of analytic functions in the variable z. Given a normalized state |<fr) = Ylk=m Ck\k, m),Ck G 
C on $) m , we obtain 



m {z\<5>) = ^r^=Y -n= (55) 
VUM) t? n y/T(n -m + l)T(n + 1) 

such that the entire functions 



- Ill 



f(z; m) = m{m = y S = = (56) 

are analytic over the whole z plane. Then, from (143]) . we can write 



|$) = f d Q {z) ^= f{z-m)\z) m (57) 

Jc y/I m (2\z\) 

and express the scalar product of two states |$i) and |$2) given on Sj m by the formula 



($i|$ 2 )= / dg(z) ' ' fx{z]m)f 2 {z-m). (5* 



s 



3.3 Mean values 

Using equation ( )35|) . the mean value of a physical observable O with respect to the BGCS 
\z) m is obtained as: 



MO\z) m = (0) z>m = -LL-r J2 * n 



m —k—m 

I m (2\z\) ~ l/r(n-m + l)r(n + l) 

n,k=m 



(k,m\0\n,m). (59) 



r(fc - m + l)T(k + 1) 
Setting n = m + ^, k = m + v, respectively, (1591) can be rewritten as follows: 



(o) g , m - , . ,s 



Im(2\z\) v j^ y/TfP + 1)I> + m + l)T(v + l)T(v + m + 1) 
x(u, m|C|v, m). (60) 

The computation of the mean values of the (i — 1, 2, 3), from (]24j) and fl60|) gives 

J +°° | z |2^+m 

(/C - } - m = Z lJ^]^T(v + l)T(v + m + l) =Z > m 

and 

(fc+)z,m = Z rr ^ — j — r = Z. (62) 

We can conclude, therefore, that (K—) z ,m an d (/C+) Z)Tn are mutually conjugated. 

Since the generators /C+, /C_ are given by /C± = JC± ± z/C2, we obtain from (JBTj) and fl62l 
together the following expressions: 

(JCi)z,m = + fc+) z > m = \( z + z) = Rez ' (63) 



1 1 

(K, 2 )z,m = ~~ fc+)z,m = ~{Z ~Z) = -ImZ. (64) 

In order to compute the mean values of the generator /C 3 and its second power /C|, it is 
useful to evaluate the sum S n , with n — 0, 1, 2, . . ., given by (see [21] in Appendix): 

= V ^ 7 (65) 

^ r (1/ + 1 r {y + n + 1 



Then, from f l24|) . we obtain 



/r v 1 | W(2M) m+1 

(/C3) - =N X(2Nr + ^ (66) 



and 



-2\ I ,2 J m+2(2k|) , , M_,-Wl( 2 kl) , fm+1 



Besides, the number operator N diagonalizing the basis vectors {\v,m),v > 0} for the 
number states, 

N\v, m) — v\v y m), (68) 
is used to compute the photon number distribution as follows: 

\(u,m\z) m \ 2 = — , Typw . —tv- (69) 

l m (2\z\) T(v + l)Y{v + 771 + 1) 

Exploiting (1661) and ( |67|) . we are able to compute explicitly the mean values of the number 
operator and its second power to obtain: 

(N) zm = (lC 3 -^±±) =\zM 2 } Z l (70) 
S \ 2 / 1 1 I m (2\z\) 1 ' 



K%-{m + 1)/C a 



m + 1 



1 2 ^m+2 



(2|*|) 



im(2|z|) 



+ z 



<m+l 



I m (2\z\) 



(71) 



Then, it becomes straightforward to extend the calculation to the intensity correlation de- 
fined as in 



( 2) = (iV 2 ) z , m - (N) Zim = I m (2\z\)I m+2 (2\z\) 
9z ' m (N)l m [W2M)] 2 



(72) 



More specifically, for two interesting limiting cases of the \z\ variable, i.e. for \z\ <C 1 and 
\z\ ^> 1, using the approximations for the Bessel modified function I m (x) [22] 

1 fX 



Im(x) 



T m + 1) V2 



respectively 



I m (x) 



l + OI - 

x 



V27TX 

we get for the intensity correlation function, with 

|_|J71 I ,,1771+1 



I m (2\z\) 



T(m + 1)' 
the following expressions: 



Im+l(2\z\) 



g(2) „ 

j z.m 



T(m + 2) 



m + 1 

m + 2' 



'm+2 



(2kl) 



I m+2 



T(m + 3)' 



(73) 



(74) 



(75) 



(76) 



g£L * i, (77) 

respectively. Thus, for small values of the intensity correlation function is smaller than 
unity, for all values of m. The corresponding BGCS \z) m exhibit sub-Poissonian statistics 
behaviour, while for large \z\, these states tend to have Poissonian statistics. In addition, 
the photon-number distribution fl69l) is sub-Poissonian [23] . 
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3.4 Probability density and time evolution 

From the quantity 



M'\ z )i 



I I \ m/2 

\Z Z\ N 



E 



(Vz'z) 



z'z J ^/I m {2\z'\)I m {2\z\) ^ I> + l)T(v + m + 1) ' 
given a normalized state |zo)m, let us define the probability density as follows: 



I in 



I m (2\z\)I m (2\zo\) 
The time evolution behavior of Q Zo (z) is then provided by 



(7S 



(79) 



By acting the evolution operator U(t) = e~i Ht on the state \z ) m , we get 
\ z o',t) m = e %Ht \zo) m 



e 2 



[n(n+oj c )+Cl]t 



m/2 



E 



e -§[n(n+o; c )+f7]t 



^IJ2\zJ) ^ v/r(z/ + l)r(z/ + m + l) 

z o(*))m 



1 1/, m) 



where Zo(t) := 2 o e 



-i[n-w c ]* 



Then, we obtain 



50) 



(81) 



feo(M) 



J m (2 v ^(t)i)/ m (2 v ^(t)i) 



7 m (2|z|)7 m (2|zb(t)|) 
It comes that the time dependence of a given BGCS \z) m is provided by 



52) 



\z;t) m = e * m \z) 



e 2 



|[n(n+w c )+n]t 



z(t)) m , z(t) := ze-^^K 



4 Quantization with the su(l, 1) coherent states 



(83) 



As is proved in Section 2, the su(l, 1) CS family resolves the unity. As an immediate 
consequence, we establish in this section the correspondence (quantization) between classical 
and quantum quantities. For more details in the quantization procedure see [2U [2S] and 
references listed therein. 

4.1 Quantization of elementary classical observables 

The Berezin-Klauder-Toeplitz quantization of elementary classical variables z and z is real- 
ized via the maps z i — > A z and z i — > A s defined by 
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A z := / z\z) mm (z\dg(z), A z := / z \z) m m (z\dg(z) (84) 



in the Hilbert space, which give, by using the equations (124| . (JbTj) and (j62|) . the following 
relations: 



A z = (fc-} Zjm \z) mm {z\dQ(z) = ^fc-\n+l,m){n + l,m\ 

^ C n=m 
+00 

= ^/(w — m + l)(n + l)|n, m)(n + 1, m\ (85) 



and 



A s = / {K + ) Zt m\z) mm {z\dg(z) = ^2 K>+\ n ~ l,m){n- l,m\ 

^ C n=m 
+00 

= yj n(n — m) |n, m) (n — 1, m| (86) 



with the matrix elements: 



OO /"27T 

/ — t[(k— to)— (n+1— m 



(A z ) n , fe = - / K m (2p)p n - m+k+2 dp x / e -#- m )-C+i-»»^ x 

Jo Jo 

x - (87) 
v/r(n - m + l)r(n + l)r(fc - m + l)r(fc + 1) 



r*oo /*2tt 

(A z ) n , k = 3 / K m (2p)p n - m+k+2 dpx / e - i[(fc+1 - m) - (n - mW rf0x 
'0 



v/r(n - m + l)T(n + l)r(fc -m + l)T(k + 1) 



respectively. 

Their commutator [A Z ,A Z ] well gives the expected result 



+00 



[A g ,A g ] = ^22K 3 \n,m){n,m\, (89) 

n=m 

keeping in mind the su(l, 1) commutation rules ( 1211) . 

Consider now the usual phase space conjugate coordinates (q,p) through z = Then, 
we get for the classical position and momentum functions q and p the Hilbert space operators: 

Q = A q :=^=(A z + A- z ), (90) 

and 
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P = A p := -±=(A Z - A- z ) (91) 



leading to a commutation rule reflecting the su(l, 1) algebra, given by 



+00 



[Q,P] = i ^2)C 3 \n,m){n,m\. (92) 

n=m 

Furthemore, using (|85|) and (|86l) . we get 



- +00 

^i (2+s) = 2E v / (^- m + 1 )( ra + 1 )l ri > m )(^ + 1 ) m l 

+00 

+ ^/ n(n — m)|n, m)(n — 1, m|], 



(93) 



and 



. +00 



= V / («-^ + l)(w+l)|n,m)(n + l,m| 



+00 



by: 



— a/ n{n — m)\n, m){n — 1, m\). 

n=m 

(94) 

Other interesting results emerging from this context are the following mean values given 

m (z\A 2 z \z) m = z 2 , m (z\A 2 z \z) m = z 2 , (96) 

m{ z \AzA z \z) m = \z\ , m(z 

\A Z A- Z \z) m = |z| 2 + 2(/C 3 ) 2 , m . (97) 

4.2 Other relevant statistics and quantum optics properties 

Now, we exploit the results issued from the quantization procedure performed above to derive 
quadrature operator mean values and some quantum optical properties of the constructed 
BGCS \z) m . 
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4.2.1 Mean values 

The mean values of the operators Q, P, Q 2 and P 2 are derived in the BGCS \z) m usin^ 

dHSD-dniD as 



l (z\Q\z) m = (Q) z , m = q, m(z\P\z)m = (P) z ,m=P, (9£ 



m(z\Q 2 \z) m = (Q 2 ) z , m = q 2 + (1C 3 ) Z;m , m (z\P 2 \ Z ) m = (P 2 ) z>m = p 2 + (K 3 ) z>m . (99) 



As well stated indeed in [24], these formulas give coherent states a quite classical face, al- 
though they are rigorously quantal! The physical meaning of the variable x in the Schrodinger 
position representation of the wave is that of a sharp position. Unlike x in the Schrodinger 
representation, the variables p and q represent mean values in the coherent states. For this 
reason both values can be specified simultaneously at the same time, something that could 
not be done if instead they both had represented sharp eigenvalues. 
There follow the dispersions: 

• For | z | <C 1 by using (|66|) and the approximation (!75|) for the Bessel modified function 



(AQ)? m = (Q 2 ) z , m - ((Q} z , m ) 2 = (/C 3 ),, m ~ 1 ' + ,T ' (100) 



2(m + 1 



yielding 



(AQ) 2 (AP) 



Z ,m 



2\z\ 2 + (m + l) 
2(m + 1) 



(102) 



which reduces for both conditions \z\ — > and m = to (AQ) 2 m (AP) 2 m ~ ~, i.e. 
the saturation of the Heisenberg inequality. This is one of the most important features 
exhibited by coherent states: their quantal face is the closest possible to its classical 
counterpart. 

For \z\ ^> 1 by using ( |66|) and the approximation ( 1741) for the Bessel modified function 



777 — I— 1 

(AQ) 2 = (Q 2 ),, m - ((Q) z , m ) 2 = (JC 3 ) z , m * \z\ + — ±— (103) 



777 — I— 1 

(AP)J ra = (P 2 ) z , m - ((P) z , m ) 2 = ()C 3 ) z , m ~ |*| + (104) 
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providing 



(AQ)'(AP) 



2 _ 



m + 1 

Z + — : 



(105) 



It is worth noticing that, in all cases, (AQ) 2 zm = (AP) 2 zm as in standard CS. This is a 
universal property exhibited by coherent states. 

4.2.2 Signal-to-quantum-noise ratio (SNR) 

For a normalized state in terms of the self-adjoint quadrature operator Q, the SNR is 
defined as [26] 



°m = jajjjr. (106) 

A low SNR indicates that the noise dominates the measurements in the state \(f>). In opposite, 
a high SNR implies that the measurements are relatively clean in the state \<p). Then, for 
q = y/2\z\ cos0, the SNR in the BGCS \z) m is given as follows: 

• For \z\ <C 1 by using (llOOp . 



((Q) z , m ) 2 q 2 
az ' m ■ (AQ) 2 (AQ) 2 



z.rn 



4:\z\ 2 COS 2 <JJ 

m + 1 , ,' ' 107 

J 2 z 2 + m + 1 2 



with u z , m — > as well as \z\ — )■ indicating an absence of noise in the measurements; 
For \z\ 3> 1 by using f |103j) . 



((Q),, m ) 2 _ ? 



(AQ)l m (AQ)l m 

4kl 2 cos 2 ^ 



2|z| + m + 1 

attesting the presence of noise in the measurements for cos0 ^ 0. 



4.2.3 Mandel parameter 

Several parameters can be introduced to characterize the statistical properties, and the most 
popular one is the Mandel parameter [28l [3Tj . denoted here by Q, known as a convenient 
noise-indicator of a non-classical field, which is frequently used to measure the deviation from 
Poisson distribution, and thus to distinguish quantum process from classical ones |28]-|35j. 
The Mandel parameter Q is defined as [28] 
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(AAQ 2 - (AQ 
(AO 

= ^ / dt 2 [ 2 dtx[l + A(tx)] - (J)T, (109) 

where (A 7 ") is the average counting number, (AN) 2 is the corresponding square variance, 
(I) = (N)/T is the steady-state photon-counting rate expressed in units of cps, and A(r) = 
(AI(t)AI(t + t))/ (I(t))(I(t + t)) is the normalized two-time correlation of intensity fluctu- 
ations (AI(t) = I (t) — {I if))) of time difference equal to r [27J . 
Moreover, the Mandel parameter Q 

2 = ^^-1 = ^-1 (110) 

is closely related to the normalized variance, also called the quantum Fano factor J 7 [32] , 
given by J 7 = (AA r ) 2 /(A r ), of the photon distribution. For J 7 < 1(Q < 0), the emitted 
light is referred to as sub-Poissonian; J 7 = 1, Q = corresponds to the Poisson distribution, 
whereas for J 7 > 1, (Q > 0) the light is called super-Poissonian [2B]-[M]- Using (ITUj) and 
11]), we have in the BGCS \z) m : 

For \z\ < 1 



z| 2 \z\ l0 kl 2 



(m + 2)(m + l) 1 'm + l 1 1 (m + l) 5 



(A^) 2m ~ \z\- u — (111' 



providing J 7 < 1 and 
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(m + l)(m + 2) 



2 T7 7 <0 (H2) 



indicating that the BGCS \z) m have sub-Poissonian statistics as discussed in the sub- 
section 13. 3( 

For \z\ > 1 



(AN)l m ^\z\, (N),, m ~\z\ (113) 

giving J 7 ~ 1 such that 



Q~0 (114) 

which implies that the BGCS \z) m have Poissonian statistics for \z\ ^> 1 as also noticed 
in the subsection I 
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4.3 The quantum energy operator 

We now deal with the energy operator for the one- dimensional quantum harmonic oscillator. 
In this case, the classical Hamiltonian H = (q 2 + p 2 )/2 = \z\ 2 should be mapped to the 
Hilbert space operator 

„ +00 
A\ z \2 := / \z\ 2 \ z) mm (z\dg(z) = ^(n - m+ l)(n+ l)\n,m)(n,m\ 



n=m 
+00 



A x A g = A g A„+ ^ 21C 3 \n,m){n,m\ (115) 



which, using (190]) and (151]) . is re-expressed as 



^ 2 + p2 , V^r- 1 , 
^M 2 = 9 h ^3|w,m)(n,m| 



2 



-00 



(^ 2 -)- .P 2 1 x 

2 (2n — m + 1) |n, m) (n, m| 



2 2 

)2 , p2 



O 2 + P 2 



It is worth noticing that, using the relations q 2 = \z\ 2 + £ -j £ ^ and p 2 = \z\ 2 — ^t** , ( 11 15H 
and the quantized values of z 2 and z 2 : 



+00 



A z 2 = ^ V( n -m + 2){n + 2)(n - m + + l)|n,m)(n + 2,m\, (117) 

n=m 

+00 

A S 2 = ^ y/(n-m + 2)(n + 2)(n-m + l)(n + l)\n + 2,m)(n,m\, (118) 

n=m 

we obtain 

+00 

A q 2 = y^(n — m + l)(n + l)|w, m)(w, m| + 

+ - x a/(ti - m + 2)(n + 2)(n -m + l)(n + l)x 

n=m 

{\n, m) (n + 2, m\ + \n + 2, m) (n, m\}] 

+00 

= Q 2 + ^/C 3 |n,m)(n,m| + |2,0)(0,0| + v / 3|3,l)(l,l| (119) 

n=m 

and 

+00 

A p 2 = y~^(n — m + l)(n + l)|n, m){n, m\ — 
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J +00 

- x Vi 71 -m + 2)(n + 2)(n -m + l)(n + l)x 



2 

n=m 

{\n,m)(n + 2,m\ + |n + 2, m)(n, m\}] 



+00 



= P 2 + ^/C 3 |n,m)(n,m| - |2,0)(0,0| - ^3, 1)(1,1|. (120) 

n=m 

From fU19p and (11201) . comparing the operators Q 2 , A q 2 and P 2 , A p 2, respectively, one remarks 
that they differ by a multiple of the commutator [A^A^] given in ( 189]) plus orthogonal 
projectors. Besides, summing the expressions f 1 1 1 9 j) and (11201) yields the relation f II 1 6 j) 
determining the operator A\ z \2. Finally, the commutators between A\m and A z , A z , A z 2, A z 2 
are calculated as follows: 



+00 

[A| 2 |2,A 2 ] = - y^(2ra - m + 3)^(7^ - m + l)(n + Tj\n,m)(n + l,m|; (121) 



[A| 2 |2,t4 5 ] = (2n — m + 1) y/ n(n — m)\n,m){n — l,m\ 

n=m 
+00 

= 2 >Jn(n — m)lC^\n,m){n — l,m|; (122) 



+00 



[v4| 2 |2,v4 z2 ] = -2 ^[2ra - m + 4]-y/(rt - m + 2) (re + 2) (re - m + l)(re + 1) 

x|n,m)(re + 2,m|; (123) 



[v4| 2 |2,v4 5 2] = 2 ^[2n - m + 4]y/(w - m + 2) (re + 2)(re - m + l)(re + 1) 

n=m 

x\n + 2,m)(n,m\. (124) 

5 Thermal properties of the BGCS 

This section furnishes a description of the statistical properties of the BGCS for the studied 
model. Consider a quantum gas of the system in the thermodynamic equilibrium with a 
reservoir (thermostat) at temperature T, which satisfies the quantum canonical distribution. 
The corresponding normalized density operator for a fixed number rre is given by 

1 00 f 1 \ H 

Pm = -J ^2e- pE ^"\u,m)(u,m\, E UiV = Ml ( n + - J - - (fi - u c ) (n - u), (125) 

where the partition function Z is taken as the normalization constant. 
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The diagonal elements of the density operator in the BGCS representation are given by 



\2u 



>,{ z \Pm \ z )m — 7 7 /ni |\ / , C ^ TV i i \ TV i i 1 \ ' (126) 

Z I m (2\z\) ^ T(v + l)T(u + m + 1) 



Note that the quantity m { z \Pm\ z )m is analog to the "semi-classical "phase space distri- 
bution function p(x,p) = (z\p\z) associated to the density matrix p (here p m ) of the system 
which is normalized as J(dxdp/27rh)fi(x,p) = 1, and often referred to as the Husimi distri- 
bution 



Then, from (1331) and (I125P together with E n ^ v = |[n(f2 + cj c ) + |(f2 — cj c )i/, we deduce 



m{z\ Pm \z) m = l-fKW] e f m( a- Uc ) ^^'^"mT^ . (127) 
The normalization of the density operator leads to 

Trp m = / dg(z) m (z\p\z) m = 1. (128) 



By the use of the following integral relation 



2 A+1 I> + l)a~ x+v+1 

1 1 x 1 1111 1 b 2 ' 
AH — u H — t;, A u H — v: v + 1; — 

2 2 2^ 2 ' 2 2 2^ 2 ' 'a 2 

+ 1 - A ± ^) >0,a>6], (129) 



we obtain 

z = e -f[n(u + ^)+a] F ^ m + i ; i ;m + i ; e -^Mj , (130) 

where F corresponds to the hypergeometric function satisfying the following property [38] : 

F(c, /i; c; x) = F(/i, c; c; x) = (1 - x)~ M . (131) 
It comes that the partition function Z takes the form: 

Z = e -f[(-+IX n +^)]|2sinh|^[Q-a; c ]}| * (132) 
and the diagonal elements of the density matrix (11271) can be written as 

MpM m = 2ef I»-l<"-'> S inh{^[ fi -.j} '"'^"J"""''' ■ (133) 



This above expression becomes, in the case of a strong magnetic field, i.e. Uq <C uj c 

lj2\z\e-*2 



m (z\ Pm \z) m = e h ^ m -V h(3^ \ ... 7 ■ (134) 

Equations (11331) and (11341) suggest that we look for another thermodynamical aspect of 
the studied model in the BGCS representation. CS are relevant to the concept of Wehrl 
entropy W gQl [391 HI] defined as 



2 'm 
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W :-- 



dxdp 
2nh 



fj>(x,p) ha./jb(x,p) 



(135) 



where fj,(x,p) = (z\p\z). Evaluating the Wehrl entropy in terms of the distribution 
function m (z\p m \z) m in (I133p . we obtain 



W 



dg(z) m (z\p m \z) m \n[ m (z\p m \z) 



-ln<{2sinh 



2e f[^c](m-i) sinh< )^^ 



4 

ijj, 



[Q — uj c ](m — 1) 



x 



f d*z 2 -K m {2\z\)U2\Ae-^^ ] ) 



(136) 



Assuming that the argument of the logarithm function under the integral is dominated 
by the region \z\ <C 1, we get approximately 



d 2 z-K m (2\z\)I m (2\z\e 

7T 



ph e ^[n- Wo] (i- TO ) 

~ — m— lit — ujA x — 3z -. 

4 L cJ 2sinh{f [O-wJ} 

Therefore, the related Wehrl entropy is approximated to 



1 { UM) J 



(137) 



W 



In 



1 - e 



-pn- 



-l + W A 



HO 



'1381 



where Who = 1 — In 



1 - e 



-pn- 



is the conventional harmonic oscillator Wehrl entropy 



[39] with frequency fi_ := „^ c . Besides, setting 



2-71-/2 

(^]p|^) m = ( — ) m {z\p\z), 



(139) 



where Z_ :- 



M- 



and A = irTZ 2 with 1Z the radius of the cylindrical body considered 



in [2], and evaluating 



W 



A 



2nl 2 



dg(z) m (z\p\z) m \n[ m (z\p\z) m ] 



(140) 



as in ([ESD-tPBD, we get 



W~-hx 



1 - e-? h -z 



, 2vr/ 2 , 

( ; 



I4r 
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where W ca i c = 1 — In [l — e @ m ~\ — In ^^j- J is the Wehrl entropy calculated for Landau's 
diamagnetism for a spinless electron in a uniform magnetic field where the thermal harmonic 



oscillator [12] frequency is denoted by f2 and / = y j^. Hence, the present physical model is 
an approximation of the problem of a thermal harmonic oscillator with frequency f2_ :— Q ~ Uc 
In the case of a strong magnetic field, from (11341) we get 



ncr'illatnr with fr^nn^nr^^r O ~2 



w 



/^-In^MY 
2uj c \ u c J 



2 sinh 



2u) c 



(142) 



The diagonal expansion of the density operator can be performed for this physical model, 
as investigated in [23] in the constructed BGCS as follows: 

p m = dg(z)P m (z)\u,m)(u,m\, (143) 
Jc 

where the function P m (z) must be determined. To this end, we determine at first the diagonal 
elements of the density operator p m in the basis of the number states, namely {\u, m), v > 0}, 
by setting: 

(v,m\p m \v,m) = / 2 K m (2\z\)I m (2\z\)P m (z){v,m\z) mm {z\v,m). (144) 

Jc 71 

By the use of J2T=o m )( v i rri| = 1 and (I125p . they are given by 

(v,m\p m \v,m) = -e~^ = [1 - e~f («-^)] e -f (n-u„> (145) 

z 

By use of the integral (|4T|) , we obtain for the function P m (z) the following expression: 

ft.M = [.^ - ll^' -'- Jf - (2 ' Z| fT )) (146) 
which is normalized as 

d e (z)P m (z) = 1. (147) 
Then, the diagonal representation of the normalized density operator takes the form 

Pm = [e^^ - l]e^^ [ d 6 (z) Km{2 ^ ^ Wc)) \z) mm {z\- (148) 

JC K mK 2 '\z\) 

Therefore, given an observable O, one obtains its mean value, i.e., its thermal average 
riven by 

(0) m = Tr{ Pm O) = J dg(z)P m (z) m {z\0\z) m . (149) 



c 



In this manner, the thermal expectation value of the number operator N, by using (1701) . 
ffT30l) and ffTBTD together, is 



1 



(N) m = Tr(p m N) = I d Q {z)P m {z) m (z\N\z) m = -. (150) 
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For the operators Q and P, setting q = y/2\z\ cos (ft and p = y/2\z\ sin0, we obtain from f|98|) 



(Q) m = 0, (P) m = 0. (151) 

In the same vein, from ( 17T1) . ( 1130f) and ( 11311) together, the thermal expectation value of 
the square of the number operator becomes 

(N 2 ) m = Tr(p m N 2 ) = [ dg(z)P m (z) m (z\N 2 \z) m : 



(152) 

One remarks that both thermal expectation values (N) m and (N 2 ) m are independent of the 
Bargmann index given here by m. One can therefore define the thermal intensity correlation 
function, which is also independent of the index m, by 

te>- - = ^ = * (153, 

In the case of Q 2 and P 2 , by use of (1991) . we have 

(Q2) m = [ e f («-«c) _ l] e -WM x F (m + 2, 2; m + 1; e"^^) x (m + 1) 

/ & h ro ,1 v 1 777- + 1 

+ ( e ^-[n-«c] _ !)-! + 



2 ' 

(P 2 ) m = [ e ^(n-«c) _ 1]e -^^ c ) x F + 2, 2; m + 1; e"^ ^) x (m + 1) 

+(e f m-«j _ i)-i + r^+i. (154) 
6 Concluding remarks 

In this work, we have investigated the Fock-Darwin Hamiltonian describing a spinless elec- 
tron gas, subject to a perpendicular magnetic field B and confined in a harmonic potential 
^Mu 2 r 2 . Thanks to the set of introduced step and orbit-center coordinate operators, we 
have shown that the studied system possesses su(l, 1) Lie algebra. As a consequence, CS 
have been constructed as the eigenstates of the 57/(1,1) group generator /C_. The mean 
values of SU(1, 1) group generators, the probability density and the time dependence of the 
BGCS have been discussed. By using these CS, the Berezin - Klauder - Toeplitz quantiza- 
tion has been performed. In this framework, mean values for quadrature operators have been 
computed highlighting the similitude between the BGCS quantal behavior and its classical 
counterpart. Furthermore, the relevance of this quantization results in deriving some impor- 
tant quantum optical properties such as the SNR and the Mandel parameter has been put 
into evidence. Statistical properties of an electron gas in the thermodynamic equilibrium 
with a reservoir (thermostat) at temperature T, satisfying the quantum canonical distribu- 
tion, have been investigated and discussed. Finally, with the density matrix provided in the 
BGCS representation, we have noticed, via the calculation of the Wehrl entropy, that the 
studied system can be identified to a model of a thermal harmonic oscillator. 
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